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We characterize vesicle adhesion onto homogeneous substrates by means of a perturbative ex- 
pansion around the infinite adhesion limit, where curvature elasticity effects are absent. At first 
order in curvature elasticity, we determine analytically various global physical quantities associated 
with adhering vesicles: height, adhesion radius, etc. Our results are valid for adhesion energies 
above a certain threshold, that we determine numerically. We discuss the haptotactic force acting 
on a vesicle in the limit of weak adhesion gradients. We also propose novel methods for measuring 
adhesion energies and we suggest a possible way of determining the size of suboptical vesicles using 
controlled adhesion gradients. 

PACS numbers: 87.16.Dg, 68.35.Np, 68.03.Cd 



I. INTRODUCTION 

When phospholipids are dissolved in an aqueous solu- 
tion, almost all the molecules condensate into bilayers. 
Lipid bilayers are formed by two contacting monolayers 
of opposite orientation, in which the hydrophilic heads 
of the molecules are located at the sides of the structure, 
the hydrophobic tails being shielded from contact with 
water 0. As there is a prohibitive energy cost associ- 
ated with their free borders, these bilayers form closed 
objects, which are called vesicles. For some biological 
studies, vesicles are used as models of the membrane of 
living cells . They also have applications as encapsula- 
tion vectors for drug delivery [0] . Their efficiency as drug 
delivery vectors is linked to their permeability, which can 
be affected by adhesion phenomena M . Vesicle adhesion 
on a solid substrate, followed by its rupture and fusion, 
also provides a simple technique for obtaining supported 
membranes Q that can be used for the design of biosen- 
sors Q. 

Adhesion phenomena between a lipid bilayer and a sub- 
strate can be divided into two categories: i) specific adhe- 
sion between aparticular host protein and a receptor on 
the substrate J7J ; this kind of adhesion generally implies a 
process of molecular recognition between a receptor and 
a ligand, and is common in biological systems, ii) non- 
specific adhesion between the membrane's lipids and the 
substrate, mediated by universal interactions, e.g., van 
der Waals forces. Here, we focus on non-specific adhe- 
sion, which can be described by an adhesion potential 
W that represents the free energy gain per unit area of 
contact. Typical values of W range from 10 _4 mJ/m 2 to 
lmJ/m 2 ||. Note that the description of adhesion us- 
ing a contact potential is approximate, because van der 
Waals forces are actually long-ranged and because mem- 
branes may fluctuate in the vicinity of the substrate: ad- 
hering vesicles actually never strictly come into contact 



with their substrate. Membrane-substrate separations 
range from 1 nm for the strongest values of W |3| , to 
about 50 nm for the weakest adhesions The highest 
values of W tend to produce vesicle rupture during the 



adhesion process 
in the membrane 



owing to a strong tension induced 



To determine the shape and free energy of an adher- 
ing vesicle, one must take into account the competition 
between: i) the adhesion energy gain, ii) the constraints 
on the total membrane area A and the total enclosed 
volume V , and iii) the free energy cost associated with 
the curvature elasticity of the membrane. The latter 
is described by a free energy density proportional to 
the square of the local mean curvature |u|. For lipids, 
the corresponding bending rigidity k is of the order of 
fO -19 J ~ 25 fee? 1 at room temperature |11 . Refined 
vesicle models take into account a constraint on the dif- 
ference between the areas of the two monolayers |f2[ |, or 
an elasticity associated with it jp| . Physically, this arises 
from the fact that lipids are not significantly exchanged 
between the two monolayers during typical experimental 
times. It is not known at the present time whether this 
constraint is significant for adhering vesicles: to simplify, 
we shall disregard it in our approach. 

The shapes of axisymmetric adhering vesicles can be 
determined by functional minimization ^ However, 
due to non-linearities in the equilibrium equations, ex- 
act solutions can only be determined numerically. In the 
asymptotic case of infinitely strong adhesion, W — > oo 
(or equivalently n — > 0), the problem is easily solved 
analytically |14j| : the equilibrium shapes are spherical 
caps, whose features are dictated by the geometrical con- 
straints only. In this paper we characterize the adhesion 
of vesicles in the case of strong but finite adhesion, by ex- 
tracting analytical corrections with respect to the infinite 
adhesion case. We determine analytically the first-order 
corrections to various physical observables and we dis- 



2 



! z 



/T x 








-7-7-7-7-7-7-7-Z-r-rv — 


r 




L li 










\ 6 







FIG. 1: Definition of the global observables associated with 
an axisymmetric adhering vesicle. The vesicle's shape, which 
was calculated numerically, corresponds to a rather deflated 
situation in which the "contact angle region" is broad and not 
well-defined. When adhesion is stronger, the vesicle's shape 
resembles a spherical cap (dashed line) , with a strongly curved 
rim at the foot of the "contact angle" (no discontinuity of the 
membrane's normal). 



and their adhesion onto homogeneous substrates. We 
also define various global observables relevant to the ad- 
hesion geometry. Section ITl| contains the results of our 



analytical calculations: in Sec. [II A , we reca ll the asymp- 
totic limit of infinite adhesion; in Sec. Ill B, we recall the 
general equations describing the equilibrium shapes of 
adhering axisymmetric vesicles; in Sec. Ill C 1 we calcu- 
late the shape of the contact angle region at first-order 
in y/n/ (WA); in Sec. Ill C 2|, w e determine the con tact 
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angle extrapolation length |17|, |20fl ; in Sec. 
determine the first-order expansions, in power series of 
y/K/(WA), of the various global observables associated 
with the vesicle's shape; in Sec. |III C 4 , we calculate at 
first-order the free energy of adhering vesicles and we 
discuss haptotaxis (motion induced by an adhesion gra- 
dient) pl[ |. In Sec. [V, we determine numerically the 
global observables and we discuss the range of validity 
of the corresponding first-order expansions. Finally, in 
Sec. [v|, we summarize our results and we discuss some 
possible applications, including new methods for mea- 
suring W. 



II. DESCRIPTION OF ADHERING VESICLES 



cuss their limit of validity by a direct comparison with 
exact numerical results. 

The first-order corrections with respect to the infinite 
adhesion limit originate from the existence of a strongly 
curved region at the border of the adhesion disc (see 
Fig. . We shall refer to this region as the "contact angle 
region" , by analogy with wetting phenomena |15fl . The 
shape of this region has been determined in Refe. Jl6|, \l% 
using an open membrane description, i.e., no volume con- 
straint and an externally imposed tension acting along a 
fixed direction mimicking the asymptotic contact angle. 
Imposing explicitly the volume constraint, we recover the 
same shape for the contact angle region. The novelty of 
our approach resides in the analytical description of the 
various observables associated with the adhering vesicle 
(height, radius of adhesion, . . . ). 

Standard measurements of adhesion potentials W are 
based on the determination of the shape of the contact 
angle region, e.g., by Reflection Interference Contact Mi- 
croscopy (RICM). Indeed, the radial curvature c of a 
detaching membrane yields W through the equilibrium 
relation c = y / 2W/n Q In practice, it is diffi- 

cult to precisely measure c, and it is more efficient to fit 
the contact angle region using RICM 0, |0). Available 
models rest however on linearized equations for contact 
angles close to ir fl7| . Our non- linear analysis allows not 
only to fit the contact-angle region and determine con- 
tact potentials even for contact angles far from tt, but 
also provides new means of determining W by measuring 
the various global observables. 

Our paper is organized as follows: In Sec. |l], we intro- 
duce the model used to describe the elasticity of vesicles 



In most experimental situations, although vesicles are 
slightly permeable to water, their volume V is strongly 
fixed by the osmotic pressure of the various solubilized 
ions to which the membrane is impermeable Q . We shall 
suppose that this volume constraint remains satisfied for 
adhering vesicles. The area A of vesicles is also fixed to 
a high accuracy: solubilized lipids are almost inexistent 
and the area-stretching modulus k s , which is of the order 
of 100mJ/m 2 ^S> W, cannot significantly affect the area 
constraint |8| . It is traditional to introduce a dimension- 
less parameter v, the reduced volume, defined by 



V 



ln(A/4nf 2 ' 



(1) 



This quantity < v < 1 describes how much the vesicle 
is deflated with respect to a sphere (v = 1). Due to the 
constraints, it is fixed. 

V and A being fixed, the free energy of an adhering 
vesicle is given by 



F 



-WA 



adh 



dA i K (ci + c 2 ) 2 , (2) 



where c\ and c 2 are the two local principal curvatures of 
the membrane, k is the Helfrich bending constant [p!o[ , 
A&dh is the area of contact between the vesicle and the 
substrate, and W is the contact potential. As discussed 
in the Introduction, we simply model the adhesion by 
an energy proportional to the contact area. In principle, 
F should also contain a Gaussian curvature term KC1C2, 
however we discard it since its integral over the mem- 
brane is constant for a given vesicle topology, according 
to the Gauss-Bonnet theorem [(221. Therefore, there are 



only two dimensionless parameters in the problem: v and 
k/(WA). 

In the whole paper, we shall restrict ourselves to ax- 
isymmetric vesicle shapes. We define the following global 
observables (see Fig. [l ) : we call H the height of the vesi- 
cle measured on the revolution axis, and L the radius of 
the adhesion disc. The adhering area is thus A ac ih = irL 2 . 
In the regime of strong adhesion, vesicles almost take the 
shape of a spherical cap. In order to precisely define a 
"contact angle" even in the case of weaker adhesion, we 
introduce the sphere which is osculatory to the membrane 
at the point intersecting the revolution axis. We call R 
its radius and 9 the angle at which it intersects the sub- 
strate (see Fig. [j]). Finally, we define the extrapolation 
length Ai as the distance between the point where the 
vesicle detaches from the substrate and the intersection 
between the osculatory sphere and the substrate. 

We shall denote throughout by the index zero all the 
quantities referring to the limit W — + oo, where the vesi- 
cle exactly takes the shape of a spherical cap. Therefore 
H , R , 9 and L are the height, radius, contact angle, 
and adhesion radius of the corresponding spherical cap. 



III. ANALYTICAL RESULTS 

Strong adhesion corresponds to the situation where the 
adhesion energy gain is very large compared to the elastic 
energy of the vesicle. Since the energy of freely floating 
vesicles is of order k ||, even for deflated vesicles, this 
condition can be expressed as 



WA > k. 



(3) 



It corresponds, for a given vesicle, to strong enough con- 
tact potentials W, or, for a given W, to large enough 
vesicles. In this situation, elasticity can be treated as a 
first-order correction with respect to the asymptotic limit 
of infinite adhesion. We shall therefore first review the 
limit W -> oo El. 



A. Infinite contact potential W 

In this case, adhesion is the only relevant contribu- 
tion to the free energy of the system, and v is the only 
dimensionless parameter of the problem. Taking into ac- 
count the two geometrical constraints, and formally set- 
ting k = 0, the shape of the adhering vesicle is deduced 
from the minimization of the following functional: 



-WA° 



adh 



E A + P V. 



(4) 



So is the Lagrange multiplier associated with the area 
constraint and Pq is the Lagrange multiplier associated 
with the volume constraint. Equation (|j) can be rewrit- 
ten as 



(X -W)A adh + Z (A-A° adh )+P V. (5) 




FIG. 2: Definition of the parameters used in the determina- 
tion of the equilibrium shape of an adhering vesicle. 



This functional is identical to that of a liquid droplet wet- 
ting a flat substrate, with the correspondence So — W — > 
7SL-7SV, S -> 7lv and P -> -AP, in which 7 SL , 7sv, 
7lv have their usual meaning and AP is the drop's excess 
pressure Jig ] . This implies that infinitely strongly adher- 
ing vesicles and liquid droplets have the same ensemble 
of equilibrium shapes, although they are described by 
different sets of physical parameters. Consequently, the 
equilibrium shapes in the asymptotic limit W — > oo are 
spherical caps. 

The major physical difference with the case of liquid 
droplets is that the contact angles are not fixed by surface 
tensions, but rather by the geometrical constraints acting 
on the vesicles. The relation between the contact angle 9$ 
and the reduced volume v [Eq. (Q)] can easily be deduced 
from simple geometry H : 



9 cos 6*o + cos 3#o 



2(2-2 cos 6» + sin 2 9 Q 



3/2 ' 



(6) 



As for the two Lagrange multipliers, they can easily be 
found by using the analogy with wetting droplets: the 
Young relation 7lv cos 9 + 7sl = 7sv yields 



So = 



W 



1 + cos 9q 

and the Laplace law AP = 27lv /Ro yields 



P. 



2S 
Rq 



2W sin 6» 
Lq (1 + cos# ) 



(7) 



(8) 



B. The equations describing finite adhesion 

Let us now consider the case of a finite contact po- 
tential W. The equilibrium shapes are those minimiz- 
ing the sum of the bending free energy and the adhesion 
free energy, subject to the area and volume constraints. 
Considering axisymmetric shapes, we parameterize their 
contour by the tangent angle ip(s), where s € [0, s\] is the 
arc- length (see Fig. ||), such that at s = the membrane 



4 



leaves the substrate and at s = Si it attains the revo- 
lution axis. Although ip(s) alone is sufficient to describe 
the vesicle's shape, it is more convenient to also introduce 
the distance r(s) to the revolution axis f23j . In the fol- 
lowing, we shall denote by a dot derivation with respect 
to s. The two principal curvatures are ci = ip (in the 
plane of Fig. 0) and ci — (smip)/r (perpendicular to the 
plane of Fig. . Enforcing the constraints by Lagrange 
multipliers, the equilibrium shapes can be obtained by 
minimizing the following functional |2j| : 

F*[r(s),^(s), Sl ] = Trr(0) 2 (X-W) 



+ / C{r,f,i>,il),i) ds, (9a) 
Jo 



where 



2nr 



1 / • sin?/; 



p 

£ + —rsmip 



+ 2irj(s) (r — cos ip) . (9b) 
Here ip{s) an d f(s) are regular functions satisfying 



#))=0, 



and r(si) = 0, 



(10) 



while r(0) = L and Si are arbitrary. The above condi- 
tions are necessary for the vesicle's shape to be closed 
and in order to avoid discontinuities of the membrane's 
normal. The parameters £ and P are the Lagrange mul- 
tipliers associated with the area and volume constraints, 
respectively. The function j(s) is a field of Lagrange mul- 
tipliers enforcing the condition r = costp for every s: this 
allows to treat r(s) and ip(s) as independent functions in 
the first variation of F* while ensuring that r(s) and ip(s) 
effectively parameterize the same shape. 
The first variation of F* can be written as 



SF* 



ds 



dC d dC\ r , , N 
dr ds dr ) 



+ SF b \ (11a) 



with 



5Fi 



dC 

dip 



(si)Sip(s 1 



dC 

dip 



(0)^(0) 



+ ^( ai )«r(l)-^(0)«r(0) 

+ 2n(^^W) r{0)Sr(0) + C( Sl )S Sl . (lib) 

The membrane's equilibrium equations are obtained by 
setting to zero the coefficients of Sip(s) and 5r(s) in SF*: 



The constraint r = cosi/', which determines the Lagrange 
field 7(5), constitutes actually a supplementary differen- 
tial equation to be fulfilled. It is worth noticing that it 
can be obtained by varying F* with respect to 7(s), since 



dC 

97 



dC 

ds \ dj 



= 



cost/;. 



(13) 



By analogy with Lagrangian mechanics, s playing the 
role of time, there exists therefore a conserved Hamilto- 
nian H, given by 123] 



U = C-ip 



2nr 



dip 

-K ( 1p 



dL 

dr 



dL 

' dj 



sin 2 ip 



■ cos ip 



p 



-rsmip 



(14) 



For an equilibrium solution, TL does not depend on s. 

The boundary equilibrium equations are obta ined by 
setting to zero the variation 5F£ in Eq. (lib). Tak- 
ing into account Eqs. ( [To|) yields Sip(0) = 0, Sip(si) = 
— ip(s\) Ssi and Sr(si) — — r{s\)bs\. Therefore 

8F* = H(si) 5si + 2vr [(S - W) r(0) ~ 7 (0)] Sr(0). (15) 

Since Ss\ and 5r(0) are independent, we obtain 



= 0, 

{T,-W)r(0) = 7(0). 



(16a) 
(16b) 



Since Ti(s) is a constant, Eq. ( |16a| ) implies Ti(0) = H = 0. 
This yields iK-0 2 (O)+7(O)/r(O)-S = 0. Hence the above 
conditions can be rewritten as 



H = 0, 



^ 2 (o) 



IF. 



(17a) 
(17b) 



Note that Eq. (17b) is the familiar curvature bound- 
ary condition for adhering membranes and thin elastic 
plates [^4| fn|. Together with Eqs. (|To|), these equations 
form the boundary conditions of the problem. Note that 
we have 5 boundary conditions for a fourth-order system 
since s\ is also an unknown. 



C. First-order corrections to the limit W infinite 



7 sin ip Pr cos ip ip cos ip sin 2ip 
= ip — — 1 = 7TJ>-' (12a) 



KT 



= 7- -K [ Ip 



1,2 



2k 
sin 2 ip 



r 2r 2 
- E- Pr sin ip. (12b) 



In order to compute the first-order corrections to the 
limit W infinite, we shall determine the shape of the con- 
tact angle region in the case of strong adhesion. To this 
aim, we first integrate once the membrane equilibrium 
equations by replacing Eq. (|12b|) by the integral condi- 
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FIG. 3: Dimensionless universal shape of the contact angle 
region. For a given contact angle #o, the actual shape, in 
units of £ = \J k(1 + cos^oVW, is the part of the curve com- 
prised between the horizontal asymptote and the substrate, 
the latter being the tangent to the curve oriented at the angle 
6o with respect to the asymptote. 



tion H = 0: 

7 sin?/) 



7 



Pr cos tp ipcosip sin (2ip) 



Kf 

r 

cos tp 
cos ip. 



2k 'i 

P 1 / 

S + — rsin-0 — — k [ ip 



2r 2 
sin 2 ip 



(18a) 

,(18b) 
(18c) 



Shape of the contact angle region 



The equilibrium problem embodied in Eqs. (18) cannot 
be solved anal ytically. As evidenced by the boundary 
condition ( |l7b| ), the width of the contact angle region 
(see Fig. ||) is of order \f k/W ; hence the condition of 
strong adhesion can be expressed as 



«1, 



(19) 



where L = r(0) is the adhesion disc's radius. This con- 
dition refines O). We therefore start with the estimates: 



r(s) 
£ 
P 



W 

— xO(l), 

K 

Lo [l + o(l)], 
So [l + o(l)], 
Po [1 + o(l)] , 



(20a) 

(20b) 
(20c) 
(20d) 



where o(l) indicates terms that tend to zero with e and 
0(1) indi cates terms of order unity. It follows that in 
Eq. ( 18b ) all the terms in the brackets are equal to 
W x 0(1) except the last one which equals W x 0(e 2 ). 
We therefore neglect it, which amounts to neglecting the 
orthoradial principal curvature (sin - 0)/r; thus Eq. 



can be rewritten as 
Lo 



7(s) 



cos ip 



[l + o(l)] 



(21) 



Plugging this expression of j(s) into Eq. (18a) and using 
Eqs. (feOj), we obtain 



sin-0 

KCOStp 

PoLqCOSi/j ipcosil) sin 2-0 



-Kip +S H —smip 



2 k 



L 



2L% 



[l + o(l)]. (22) 



All the terms in this equation are equal to Wk~ x x 0(1), 
except the last two term which are equal to W k~ y x 0(e) 
and Wk~ x x 0(e 2 ), respectively. Using the expressions 
of the zeroth-order Lagrange multipliers ([?]) and (||), we 
obtain finally 



'. ' = ( — -ip 2 tanV> 



W sin ip — sin 9o 
k (1 + cos6> ) cosip 



o(l)]. 
(23) 



Neglecting the o(l) term provides us with a simplified 
equation describing the contact angle region in the regime 
of strong adhesion. This equation can easily be in- 
tegrated once by introducing the intermediate vari able 
tp 2 /{2cosip) and using the boundary condition (|l7b]): 



2W 1 + cos (6> + ip) 
k 1 + cos 9o 



(24) 



Its solution is 



ip(s) = 4 arctan 



tanh 



w 



4k (1 + cos6> ) 



(25) 



where we have shifted the arc-length s by a constant, the 
detachment point still corresponding to ip = 0. Since 
the radius L of the adhesion disc has disappeared, the 
problem has actually become two-dimensional, as if the 
rim of the contact angle region were translationally in- 
variant. This implies that in the present regime of strong 
but finite adhesion, the size of the vesicle has no influence 
on the shape of the contact angle region. Yet, the con- 
straint on the reduced volume keeps an influence since it 
determines #o- 

Scaling lengths to I — \J k(\ + cos6* )/VF, and intro- 
ducing a normalized frame (A, Y) rotated of an angle 9o 
with respect to the frame (r, z), the shape of the contact 
angle region assumes the universal expression: 



X(S) 
Y(S) 



2 tanh S - S, 

2 [i- (cosh sy 



(26a) 
(26b) 



where S = s/£ is the normalized arc-length. For a given 
contact angle 9q, the actual shape of the contact angle 



6 



region is obtained by putting the substrate tangent to 
this shape, at the angle 6*o with respect to the horizontal 
asymptote of the curve, and then rescaling lengths with 
respect to £, as shown in Fig. ||. This shape is the same as 
that found in which was established using an open 
membrane description and by imposing the asymptotic 
direction of the membrane at the angle #0 through an 
externally imposed tension. 



2. Contact angle extrapolation length 

A useful characteristic of the contact angle region, used 
in RICM experiments in order to determine the ratio 
W/k |t7|, ^(|, is the extrapolation length \\ (see Fig. [I]). 
From the above calculation, valid in the regime of strong 
adhesion, we deduce 



Ai 



cos ip ds 



' sin?/' 
tan 8q 



ds 



<2k a 
■ cot — . 
W 2 



(27) 



This expression holds even for deflated vesicles and agrees 
with the expression previously obtained in Rcf. 17 for 
nearly spherical vesicles (71 — 6*0 -C !)• 



3. First-order corrections to the global observables 

Let us determine, in the regime of strong but finite ad- 
hesion, the global observables characterizing the vesicle's 
shape: 0, R, L, H (see Sec. || for their definitions). To 
this purpose, we match the contact angle region to the 
rest of the vesicle. This is done by expressing the area 
and volume constraints: 



A 
V 



A 



cap 



cap 



- SA, 

sv, 



(28a) 
(28b) 



where A cap = ttR 2 [2(1 — cos 9) + sin 2 8] is the area of the 
spherical cap osculatory to the vesicle plus the area of its 
bounding disc, V C3j p = ^irR 3 [2(l — cos8) — sin 2 dcosd] is 
the volume enclosed by this spherical cap, and A and V 
are the actual vesicle's area and volume, respectively. 

In the regime of strong adhesion, SA can be evaluated 
from the results of Sec. Ill C 1 by calculating the differ- 



ence between the area associated with the approximate 
contact angle shape given by Eq. ( p5[ ) and that associated 
with its asymptote: 



SA ~ 2ttL [ ds- 2wL ( Xi + f 

Jo V Jo sm °a 



ds 



= 47 r(cos|-cot|) ^L =AxO(e),(29) 

where Lq = Ros'mOo- As for SV as SAXi, it follows 
that it is equal to V x C(e 2 ) since Ai = L x 0(e) [see 
Eq. fl27|)]. Note also that since in the limit W — > 00 
the vesicle's shape is actually a spherical cap, we have 
A = A cap {R , O ) and V = V cap (Ro, ). 




FIG. 4: Co effici ents lo (solid line) and h (dashed line) of the 
expansion (ila) of the radius L of the adhesion disc, as a 
function of the reduced volume v of the vesicle. 



Setting 9 = 8 + SO and R = R + SR, we obtain 
the first-order corrections S9 and SR by solving the sys- 
tem (p8|) to first order in e. This yields 



SR = 



2 ( sin — — 1 j (2 + cospo 

i?o sin 6*o 
2 I 1 - sin y ^ sin 2 

(1 — cos Oo) 2 




+ 0{e 2 ), (30a) 



>ZZ+O(R e 2 ). (30b) 



These results show that, in order to compensate the area 
cost SA of the contact angle region, the vesicle's shape 
flattens (SR < 0) with respect to the asymptotic case of 
infinite adhesion. 

We are now able to determine the fist-order correc- 
tions to L, the radius of the adhesion disc, and to H, 
the height of the vesicle. Since the intersection be- 
tween the substrate and the osculatory spherical cap, de- 
scribed by [R,0), is a circle of radius i?sin6*, we have 
L ~ RsinO - Ai. Using Eqs. @ and @, L can be 
written in the dimensionless form: 



L 

72 



lo + h 



WA 



0{ 



WA 1 



with, 



(31a) 



(31b) 



(31c) 



Note that 0o is linked to the prescribed reduced volume v 
of the vesicle through expression (^j|) . In Fig. [| we have 
plotted lo and l\ as a function of v. 
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h , hi 




FIG. 5: Co effic ients ho (solid line) and hi (dashed line) of the 
expansion ( B2a) of the total height of the vesicle, as a function 




of the reduced volume v of the vesicle. 



FIG. 6: Co effici ents fo (solid line) and /i (dashed line) of the 
expansion ( 35a ) of the free energy of the vesicle, as a function 
of the reduced volume v of the vesicle. 



As for H, since the osculatory spherical cap is tangent 
to the top of the vesicle, we have simply H — i?(l — cos6*). 
Using Eqs. (BO) we obtain 



with, 



ho 
hi 



' 1 — cos Oo 
7r (3 + cos 6> ) ' 



-2\/2 [ 1 - Silly 



(32a) 

(32b) 
(32c) 



The plots of ho and hi as a function of v are shown in 
Fig. g. Note that ho and lo stem from simple geometrical 
considerations, while h\ and l\ originate from curvature 
elasticity effects. 



4- Free energy of adhering vesicles 

We now turn to the determination of the analytical 
development of the total free energy of the vesicle: 



F = ~-kL z W + F el 



(33) 



where F c \ is the curvature free energy. The latter is 
the sum of a contribution F e i,i arising from the contact- 
angle region and a contribution F e i,2 arising from the 
top spherical cap. Since both the size and the cur vature 
radius of the contact-angle region are of order y/n/W, 
F e i,i is of order L^k~JW x k{^W[k) 2 = WL 2 x 0(e). 
As for -F e i,2) it can be neglected since it is of order 
k x O(l) = WL 2 x 0(e 2 ), as for a free vesicle. 

In the strong adhesion regime, the orthoradial curva- 
ture (sin?/>)/r of the contact-angle region is negligible as 



justified in Sec. Ill CI. Therefore, using Eq. (|24|), we 
obtain 



F e i ^ ttkLo / ip 2 ds = 2ttL v2kW 
Jo 



1 u 
1 — sm — 

2_ 

Oo ■ 
cos — 



(34) 



Using the expression of L given by Eq. (31a), we finally 
obtain 



F 
WA 



fa + U 



WA 



WA" 



with 



fo 
ft 



1 + COS 6q 

3 + cos 9o ' 



1 • "0 

1 — sm — 
2 

V3 + cos 



(35a) 

(35b) 
(35c) 



The plots of fo and fx in terms of the reduced volume v 
are shown in Fig. ^. 

As an application of this result, let us determine the 
force acting on an adhering vesicle in the presence of weak 
adhesion gradients: haptotaxis If the dynamical 

deformations during the movement are weak, the shape of 
the vesicle can be assimilated to its equilibrium shape on 
a substrate with a constant adhesion potential W equal 
to the average of W in the real adhesion disc. The force 
exerted on the vesicle is then 



dF 
dW 



, . 1 , I KA ( K \ 



VW, (36) 



where V is the gradient on the substrate. Since fo and 
/1 have opposite signs, the curvature elasticity decreases 
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FIG. 7: Numerically calculated radius L of the adhesion disc 
as a function of k/(W A) (solid lines) for vesicles of reduced 
volume ii = 0.77 and v — 0.93, along with its asymptotic 
expansion ( J31^ ) (dashed lines). The vertical bars indicate the 
threshold at which the relative error between the exact value 
of L and its analytical estimate reaches 5%. 



the haptotactic force with respect to the infinite adhesion 
limit. Moreover, for a given VW the haptotactic force is 
not constant but actually increases with W. 

To check the order of magnitude of the haptotactic 
force, let us consider a 10 /im vesicle (A ~ 10 _9 m 2 ) with 
k ~ 10 -19 J, subject to a contact potential varying uni- 
formly from W ~ 10~ 4 mJ/m 2 to!T~ 10 -3 mJ/m 2 on a 
distance ~ 1mm. Assuming a reduced volume v = 0.77 
corresponding to 9q ~ ir/2 [see Eq. dm], we obtain a force 
varying from 0.26 pN to 0.29 pN (8% variation). With a 
simple Stokes law, this corresponds to velocities of the 
order of 1/jms" 1 . Note that in infinite adhesion this 
gradient would give rise to a force equal to 0.3 pN. 



IV. COMPARISON WITH THE EXACT 
NUMERICAL RESULTS 

We expe ct the asymptotic expansions given in Sees. 



Ill C 3 and [II C 4 to be accurate in the regime of strong 
adhesion. To check their validity, we have compared 
them with the exact values of the vesicle's observables, 



obtained by numerically integrating Eqs. ( |l2j ) and (13). 

In order to avoid numerical instabilities when ap- 
proaching the axis of revolution (r = 0), we have chosen 
to integrate the equations starting from the top of the 
vesicle (s — si, see Fig. ||). To this aim, we impose the 4 
initial conditions: 



n{ Sl ) 

r(si) 



0. 
0. 
0. 

CO, 



(37a) 
(37b) 
(37c) 
(37d) 




FIG. 8: White area: region of the (v, n/WA) plane where the 
analytical estimate of the radius L of the adhesion disc differs 
from its exact numerical evaluate by less than 5%. Above the 
dashed line (not fully shown) the axisymmetric oblate shapes 
correspond to an unphysical self-crossing of the membrane. 



where TL{s) is the first integral of the equilibrium equa- 
tions given by Eq. <Q) and Co is an arbitrary initial 
curvature. The integration proceeds backwards, start- 
ing from s — si (the actual value of s± is arbitrary), and 
is stopped when ip — 7r, meaning that the substrate has 
been reached. To span more easily all the values of the di- 
mensionless parameter WA/k for a given reduced volume 
v, we proceed as follows. For a given fixed value of the 
initial curvature Co, we vary the Lagrange multiplier E in 
Eqs. ( |l2| ) until the solution has the desired reduced vol- 
ume. During this search, the other Lagrange multiplier, 
P, is fixed to a value (positive for weakly adhering vesi- 
cles and negative for strongly adhering vesicles) assuring 
that the size of the vesicle is of order 1 in dimensionless 
units. Once the correct value of E has been obtained, 
we determine the area of the vesicle and its curvature at 
the point ip = tt, where it touches the substrate. The 
corresponding value of W A/k is obtained through the 
boundary condition (17b). The set of all the solutions for 
a fixed v and all values of WA/k corresponds to a tra- 
jectory in the (E,co) plane that has to be reconstructed 
by varying cq and E. Sometimes, a given cq corresponds 
to two or more values of £, which yields different values 
of WA/k for the same reduced volume. 

To exemplify our results, we show in Fig. |7 the radius L 
of the adhesion disc, as a function of the reduced inverse 
adhesion energy k/(WA), for vesicles of reduced volume 
v = 0.77 (or 8q ~ tt/2) and v — 0.93. High adhesion 
energies correspond to low values of k/(WA), where our 
asymptotic formula ( |3l] ) closely fits the exact numerical 
results. The vertical bars indicate the threshold above 
which the error associated with the analytical approxi- 
mation is larger than 5%. At this threshold, L differs 
nonetheless from its infinite adhesion limit Lq by more 
than 30%: significant deviations from the infinite adhe- 
sion limit are therefore predicted with a good precision 
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FIG. 9: Same as Fig. ^ but for the total height H of the 
vesicle. 



by the asymptotic formula (pT|). 

For vesicles of reduced volume in the range 0.25 < v < 
0.95, we have determined the threshold for k/(WA) at 
which the relative error between our analytical approx- 
imations and the exact results reaches 5%. In Fig. || 
we show this threshold for the adhesion disc's radius L, 
in Fig. H for the total vesicle's height H, and finally in 
Fig. [f(J for the derivative dF/dW of the free energy with 
respect to the adhesion energy. The latter quantity is 
linked to the haptotactic force (|36|). 

Typically, the 5% threshold occurs for values of 
k/(WA) comprised between 10 -2 and 10~ 3 (see Figs. 
|l0| ). Let us consider the case of "giant vesicles" since 
they are optically observable (typical size ~ 10-100 /im). 
Supposing an area of ~ 10 3 /im 2 and a bending rigidity 
k ~ 10 -19 J, the 5% threshold occurs for values of W in 
the weak adhesion range 10 -5 -10 -4 mJ/m 2 . Our analyt- 
ical estimates seem therefore able to describe the adhe- 
sion of giant vesicles up to the lowest values of W exper- 
imentally accessible. For smaller vesicles, the threshold 
is more limitative, as it corresponds to higher adhesion 
energies W. Note also that in the case of weak adhesion, 
the picture could be quantitatively different for vesicles 
filled with a fluid denser than the outside medium, be- 
cause of gravity effects. 



V. DISCUSSION AND POSSIBLE 
APPLICATIONS 

Taking into account the effect of membrane elastic- 
ity to first-order in y/ n/(WA), we have analytically de- 
termined the global observables characterizing adhering 
vesicles. Our calculation is based on the fact that if adhe- 
sion prevails over elasticity, most of the elastic contribu- 
tions to the free energy are located in the "contact angle 
region." We have numerically determined the region of 
validity of our analytical expansions, in the (v, k/WA) 



FIG. 10: Same as Fig. | but for the derivative dF/dW of 
the free energy of the vesicle with respect to the adhesion 
energy W. 



parameter space, corresponding to a 5% maximum er- 
ror. It turns out that for "giant vesicles" (typical radius 
10-100 /im), this region comprises practically all the ac- 
cessible adhesion surface energies W . Besides, our ana- 
lytical estimates correctly describe significant deviations 
with respect to the infinite adhesion limit. 

We have throughout assumed that the area A and vol- 
ume V of the vesicle were strictly fixed, while vesicles 
actually possess small but finite stretching elasticity and 
osmotic compressibility. It is easily shown, however, that 
a self-consistent choice of the Lagrange multipliers £ and 
P yields the same equilibrium solution in the presence of 
arbitrary stretching and osmotic potentials. It follows 
that our expressions of SO, 5R, L and H remain correct 
provided that A and V are the actual area and volume 
(that now depend on W). 

Measurements of the contact potential W are usually 
performed by RICM imaging of the contact angle re- 
gion |l7], |20) . The value of W is inferred either from the 

17b]), or from the ex- 
(an approximated 



local curvature ip(0) through Eq 
trapolation length Ai through Eq. (2 
formula valid for n — 9q <C 1 is actually used |17|). The 
precision of the former measurement is limited by the fact 
that the vesicle's curvature varies abruptly close to its 
detachment point, the exact position of which is always 
slightly ambiguous. The extrapolation length measure- 
ment relies on the existence of a well defined asymptote 
of the vesicle's profile close to the contact angle region: 
it is therefore suitable only for the stron gest ad hesions. 
The expressions of L and H found in Sec. Ill C 3| allow to 
envisage novel measurements of W, based on global char- 
acteristics of the adhering vesicle. To this aim, one needs 
to know also the vesicle's total area and volume. They 
can be either directly determined by imaging a side view 
of the vesicle |24[[, or inferred by osmotically deflating a 
spherical vesicle of known radius in a controlled way. The 
fact that W can be determined through two independent 
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measurements (L and H) allows to better estimate the 
experimental errors and to validate the model. More- 
over, such global measurements are complementary to 
the above-cited local ones, since they are more adapted 
for measuring weaker values of W . The precision of the 
measurement should increase as W decreases, as long as 
one remains inside the authorized zone of Figs. |and|. 
In fact, for W too strong, L and H saturate, while for W 
too weak the analytical expansions of L and H lose their 
validity. However, as we have seen, the low W limitation 
is not relevant for "giant vesicles" . 



Finally, the haptotactic force (|36|) suggests the pos- 
sibility to determine the size of suboptical vesicles by 
measuring their velocity of migration on a substrate pre- 
senting a controlled adhesion gradient, supposing a linear 
viscous friction law. Fitting the evolution of the vesicle's 
velocity as a function of W allows to determine the vesi- 
cle's area and volume, provided that the dependence of 
the friction coefficient on n/(WA) and v is known. The 
latter could be determined using giant vesicles of known 
area and volume. 
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